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Abstract 

The hydrokinetic model is applied to restore the initial conditions and space- 
rs ' time picture of the matter evolution in central Au+Au collisions at the top 
^ . RHIC energy. The analysis is based on the detailed reproduction of the pion 
and kaon momentum spectra and femtoscopic data in whole interval of the 
transverse momenta studied by both STAR and PHENIX collaborations. A 
^ ' good description of the pion and kaon transverse momentum spectra and 
interferometry radii is reached with both initial energy density profiles moti- 
vated by the Glauber and Color Glass Condensate (CGC) models, however, 
at different energy densities. 
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1 Introduction 



In this letter we apply the HydroKinetic Model (HKM) proposed in [H [2] and 
■ further developed in [3l HI [5] to an analysis of the femtoscopic measurements 

at RHIC for central Au+Au collisions at the top energy y/s = 200 AGeV. 
Namely, we analyze pion and kaon transverse momentum spectra and the mx- 
2 I behavior of the pion and kaon interferometry radii to clarify, in particularly, 

how these observables depend on the initial conditions: Glauber and CGC- 
like. 

Hydrokinetic approach combines the advantages of the hydrodynamic 
^ ! approximation, where possible phase transitions are encoded in the corre- 

I sponding equation of state (EoS), and microscopic approach, accounting for 

a non-equilibrated process of the spectra formation due to gradual particle 
liberation. The dynamical decoupling is described by the particle escape 
probabilities in inhomogeneous hydrodynamically expanding systems in the 
way consistent with the kinetic equations in the relaxation time approxima- 
tion for emission function. The basic hydrokinetic code, proposed in is 
modified now to include decays of resonances into the expanding hadronic 
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chemically non-equilibrated system and, based on the resulting composition 
of the hadron- resonance gas at each space-time point, to calculate the equa- 
tion of state (EoS) in a vicinity of this point. The obtained local EoS allows 
one to determine the further evolution of the considered fluid elements. 

2 Model description and results 

Let us briefly describe the main features of the model. Our results are all 
related to the central rapidity slice where we use the boost-invariant Bjorken- 
like initial condition in longitudinal direction. 

Initial conditions: Two models of initial conditions for hydrodynamic 
expansion were used: Glauber model and CGC-like model. In former one, 
initial energy density in the transverse plane is proportional to the participant 
nucleon density[S]. For CGC model, the procedure described in |31[7] give us 
the energy profile in the transverse plane: 

e{xT) = eo jj^ , (1) 

Po 

where the number of participants p is defined in the same way as in Glauber 
model. The parameter eo = e{b = 0, = 0) - the maximal energy density at 
the initial moment of thermalization - is the first fitting parameter of HKM 
for both IC cases. 

Following the ideas about the pre-thermal flow development [7] we take 
a "conventional" proper-time of thermalization, Tj =1 fm/c and non-zero 
transverse flow already developed at the thermalization time. The initial 
transverse rapidity profile is supposed to be linear in radius r^: 



Vt = where Rt = J < rf, >, (2) 

here a is the second fitting parameter in our model. Note that the fitting 
parameter a should include also a positive correction for underestimated re- 
sulting transverse flow since in this work we did not account in direct way for 
the viscosity effects [6] neither at QGP stage nor at hadronic one. In formal- 
ism of HKM [2] the viscosity effects at hadronic stage are incorporated in the 
mechanisms of the back reaction of particle emission on hydrodynamic evo- 
lution which we ignore in current calculations. Since the corrections to trans- 
verse flows which depend on unknown viscosity coefficients are unknown, we 
use fitting parameter a to describe the "additional unknown portions" of 
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flows, caused both factors: by a developing of the pre-thermal flows and the 
viscosity effects in quark-gluon plasma. 

Equation of state: At high temperatures corresponding to the QGP phase 
and crossover transition to hadron phase we use a realistic EoS [11] adjusted 
to the lattice QCD results for zero barionic chemical potential so that it 
is matched with an ideal chemically equilibrated multicomponent hadron 
resonance gas at = 175 MeV. To take into account a conservation of the net 
baryon number, electric charge and strangeness in the QGP phase, we make 
corrections to thermodynamic quantities for nonzero chemical potentials, as 
proposed in [T2] . 

The chemical freeze-out temperature T^h = 165 MeV is used, with corre- 
sponding chemical potentials fiB =29 MeV, fis =7 MeV, fiE =-1 MeV and 
also the strangeness suppression factor 75 = 0.935 which are dictated by 
200A GeV RHIC particle number ratios analysis done in the statistical model 
[91 [10]. 

At the chemical freeze-out temperature T^h the EoS is matched with ideal 
Boltzmann hadronic resonance gas which includes = 359 hadron states 
made of u, d, s quarks with masses up to 2.6 GeV (the same set of hadrons 
as used in [l3]). The EoS used at temperatures below chemical freeze-out is 
chemically non-equilibrated and takes into account the change of chemical 
composition due to decays of resonances in the form: 

d^{ni{x)u>'{x)) = -TiUi^x) + ^bijTjnj{x) (3) 

j 

here one neglects a thermal motion of the resonance j, that can be justi- 
fied because post (chemical) freeze-out temperatures are much less than the 
mass of the lightest known resonance. Also, equations for the hydrodynamic 
evolution are written under supposition of an instant thermalization of the 
decay products, that is consistent with the ideal fluid approximation (mean 
free path is zero). 

Spectra formation: During the hadron phase, hadrons are permitted to 
leave the system. The process of particle emission is described by the means 
of emission function, which, for (quasi-) stable particles is expressed through 
the gain term GP™(x,p) in Boltzmann equations and escape probabilities 
[H [2]. In this formalism, particle emission is formed by the particles that 
undergo their last interaction or already free initially. Thus, emission depend 
on the collision rate of given particle with hadron medium, which is calculated 
in a standard way using the total hadron cross-sections taken in a way similar 
to that in UrQMD code. 
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Results and discussion: The results of the HKM for the pion and kaon 
spectra, interferometry radii and Rout/ R side ratio are presented in Fig. 1. 
The best fit for the Glauber IC is reached at the following values of the two 
fitting parameters related to the proper time r = 1 fm/c: eo = 16.5 GeV/fm'^ 
((e) = 11.69 GeV/fm^) and parameter of the initial velocity defined by ([2]), 
a = 0.248 {{vt) = 0.224). In the case of the CGC-like initial conditions 
Rt = 3.88 fm, the fitting parameters leading to the best data description 
are eo = 19.5 GeV/fm^ ((e) = 13.22 GeV/fm^) and a = 0.23 {{vt) = 0.208). 
Since the temperature and baryonic chemical potential at chemical freeze- 
out, which are taken from the analysis of the particle number ratios j9], are 
is more suitable for the STAR experiment, the HKM results for kaon spectra 
are good for the STAR data but not so much for the PHENIX ones. Note 
also that, in spite of other studies (e.g., jS]), we compare our results for 
the interferometry radii within the whole measured interval of pt covered 
at the top RHIC energy. Finally, one can conclude from Fig. 1 that the 
description of pion and kaon spectra and space-time scales is quite good for 
both IC, the Glauber and CGC. It is worth noting, however, that the two 
fitting parameters a and eo are various by 10-20% for different IC, as it is 
described above. 

The special attention acquires a good description of the pion and kaon 
longitudinal radii altogether with Rout/Rside ratio, practically, within the 
experimental errors. Such an achievement means that the HKM catches 
the main features of the matter evolution in A+A collisions and correctly 
reproduces the homogeneity lengths in the different parts of the system which 
are directly related to the interferometry radii at the different momenta of 
the pairs [T5] . 

3 Conclusions 

In this letter, we show how hydro-kinetic model is applied to restore the 
initial conditions and space-time picture of the matter evolution in central 
Au+Au collisions at the top RHIC energy. The analysis, which is based on a 
detailed reproduction of the pion and kaon momentum spectra and measured 
femtoscopic scales, demonstrates that basically the pictures of the matter 
evolution and particle emission are similar at both Glauber and CGC initial 
conditions (IC) with, however, the different initial maximal energy densities: 
it is about 20% more for the CGC initial conditions. The main factors, 
which allows one to describe well simultaneously the spectra and femtoscopic 
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Figure 1: The transverse momentum spectra of negative pions and negative 
kaons, all calculated in the HKM model. The comparison only with the STAR 
data are presented in the separate small plots. (Top). The interferometry 
radii and Rout/Rside ratio for ti~tt~ pairs and mixture of K~K~ and K^K^ 
pairs. (Middle and bottom). The experimental data are taken from the 
STAR [ISl [n] and PHENIX [Ml HH EO] Collaborations. 
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scales are: a relatively hard EoS (crossover transition and chemically non- 
equilibrium composition of hadronic matter), pre-thermal transverse flows 
developed to thermalization time, an account for an "additional portion" of 
the transverse flows due to the shear viscosity effect and fluctuation of initial 
conditions, a correct description of a gradual decay of the non-equilibrium 
fluid at the late stage of expansion. 
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